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Abstract —  The  limitation  of  Steam  Injection  to  depths  has 
been  a  subject  of  concern  in  the  application  of  Steam 
Injection  for  heavy  and  extra  heavy  oil  recovery.  This  is 
usually  as  a  result  of  the  complex  mechanism  of  heat  loses 
occurring  in  the  wellbore  and  consequently  the  heat  loss 
distribution  in  the  reservoir.  A  conventional  approach  to 
the  optimization  of  steam  injection  has  been  based  on 
isolated  analysis  of  the  well  system  aimed  at  maintaining 
adequate  steam  quality  at  the  sandface  at  optimal  injection 
rate,  pressure,  temperature  and  overall  heat  transfer 
coefficient.  This  often  results  to  total  neglect  of  the  effect  of 
the  interaction  between  the  well  system  and  the  reservoir 
system  in  the  Model  results.  This  research  presents  an 
integrated  approach  in  the  modelling  of  steam  injection 
project  that  incorporates  both  the  well  system  and  the 
reservoir  system.  In  this  study,  a  three  case-study  wells 
were  analyzed  which  are  located  at  INJ1  (1,1),  PR0D1  (5, 
5)  and  PR0D2  (9,  1)  respectively.  The  results  of  the 
findings  reveals  that  the  conventional  practice  of 
maintaining  sufficient  SQ  at  the  sandface  is  not  the  last 
optimization  strategy  in  real  field  scenario.  This  is  because 
the  efficiency  of  the  heavy  oil  displacement  by  the  steam  is  a 
co -function  of  the  effective  SQ  at  the  sandface,  the 
FHLRJFHLT  and  the  relative  distance  of  the  injector(s) 
from  the  producer(s)  which  are  characterized  by  the 
thermal  properties  of  the  reservoirs.  As  part  of  the 
objectives  of  this  study,  a  novel  numerical  approach  using 
PROSPER  wellbore  simulator  is  presented  for  analysing 
the  impact  of  reservoir  back  pressure  on  the  estimated  SQ. 
The  results  as  presented  in  the  work  shows  that  wrong 
estimations  of  downhole  SQ  can  result  from  the  total 
neglect  of  Reservoir  Pressure  especially  in  relatively  deeper 
wells. 
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I.  INTRODUCTION 

The  global  rise  in  oil  price  and  the  increasing  worldwide 
energy  demand  are  clear  indications  that  many  proved 
undeveloped  hydrocarbon  reserves  has  to  be  developed 
using  available  technology.  These  reserves  have  been 
identified  as  either  conventional  or  unconventional  based 
on  its  source.  The  OPEC  Annual  Statistics  (2017)  reported 
the  world’s  total  proven  reserves  as  of  2016  to  be 
1 ,492,1 64Mbbl.  According  to  PetroWiki  (2017) 
classification,  the  unconventional  sources  of  hydrocarbons 
include  heavy  oil,  extra-heavy  oil  and  Bitumen  amounting 
to  a  total  of  9  trillion  barrels  of  oil  (from  both  conventional 
and  unconventional  sources),  accounting  for  about  83.42% 
of  the  world’s  total  proven  reserves.  The  efficient  operation 
of  steam  injection  requires  the  injection  of  steam  of 
sufficient  quality  at  sufficient  rates.  However,  the  cost  of 
generating  steam  is  quite  high  making  up  about  a  half  of  the 
overall  cost  of  running  the  whole  operation.  Hence,  the 
optimized  use  of  the  injected  steam  has  been  the  industrial 
practice  (Hong,  1994).  For  optimal  application  of  steam, 
the  reservoir  depth  must  be  duly  considered  as  this  poses  a 
constraint  to  the  efficiency  of  the  operation. 

The  limitation  of  steam  injection  thermal  EOR  to  depths 
not  more  than  5,000ft  is  due  to  significant  heat  losses  in  the 
wellbore,  the  formation  and  consequently,  steam  quality 
reduction.  The  development  of  models/  simulators  is  an 
important  optimization  tool  in  a  more  modem  industrial 
society  today.  Thus,  this  enables  the  utilization  of  computer 
assisted  numerical  methods  for  the  optimization  of  the 
parameter  of  interest  over  any  possible  number  of  ranges 
for  convergence.  Hence,  an  integrated  model  that  can 
compensate  for  injection/reservoir  pressure  effect,  choice  of 
completion  design  and  the  reservoir  response  to  steam  will 
certainly  be  highly  invaluable  in  the  design  and 
optimization  of  steam  injection.  This  research  provides  an 
integrated  ECLEPSE-PROSPER  Steam  Injection  Model  for 


www.iiaers.com 


Page  |  110 


International  Journal  of  Advanced  Engineering  Research  and  Science  (IJAERS)  [Vol-5,  Issue-11,  Nov- 2018] 

httos:/ Zdx.doi.ora/1 0.221 61/iiaers.5.1 1.16  ISSN:  2349-6495(P)  /  2456-1908(0) 


steam  measurements  along  the  injection  as  well  as  the 
reservoir  response  to  steam  thermal  energy. 

Thermal  processes  are  generally  classified  as  those  EOR 
methods  that  involve  the  introduction  of  external  heat 
energy  into  the  reservoir  to  heat  up  the  high  viscous  crude 
in  the  reservoir  and  as  such  make  it  more  mobile.  The 
temperature  dependence  of  viscosity  is  of  empirical  basis. 
This  forms  the  basics  of  every  thermal  recovery  procedures 
since  the  entire  aim  is  to  raise  the  reservoir  temperature  for 
viscosity  reduction.  The  viscosity  of  liquids  as  a  function  of 
temperature  can  be  estimated  using  any  of  the  following 
correlations : 

•  The  Andrade ’s  exponential  correlation 

(Li) 

Where; 

p=  dynamic  viscousity,  cp 
T=  absolute  temperature,  K 

A  and  B  =  constants  which  varies  from  liquid  to  liquid. 

•  The  Braden  s  correlation  for  oil 

log(v2  +C)  =  (lYlogCVi  +  C)  (1.2) 

T2 

Where; 

Ti  and  T2  =  absolute  temperature  at  the  original  condition 
and  the  final  conditions  (when  temperature  is  raised) 
respectively 

Vi  and  V2=  kinematic  viscosity ,cSt 
C=  constant  (equal  to  0.6  for  io>1.5  cSt) 

D=  constant  of  the  order  3.5  to  4  (Latil,  1980) 

Generally,  hot  fluid  injection  can  be  classified  as  hot  water 
injection,  cyclic  steam  injection  (also  known  as  ‘huff  and 
puff)  and  direct  steam  injection  (also  known  as 
steamflooding),  (Latil,  1980).  For  the  scope  of  this  study, 
the  attention  is  going  to  be  concentrated  on  steam  The 
cyclic  steam  injection  also  known  as  steam  stimulation  or 
the  huff  ‘n’  puff  is  a  practice  that  uses  a  single  well 
alternately  as  injector  and  producer  for  a  more  efficient 
utilization  of  the  heat  injected.  It  basically  involves  three 
phases  of  operation  for  a  given  cycle: 

•  The  steam  injection  phase(  which  is  similar  in 
operation  with  normal  direct  steam  injection  i.e. 
steamflooding) 

•  The  soak  period  and 

•  The  production  phase. (Latil  1980) 

In  the  assessment  of  the  efficiency  of  the  steam  injection 
design,  the  major  optimization  criteria  are  to  maintain 
optimum  steam  quality  at  a  sufficient  injection  rate  using 
the  ‘ rule  of  thumb  (Hong,  1994).  However,  such  rate  must 
economically  be  considerable  to  compensate  for  the  high 
cost  of  steam  generation.  For  most  practical  consideration, 
the  ‘the  rule  of  thumb  ’  is  to  maintain  an  injection  rate  of 
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1.5B/D  cold  water  equivalent(CWE)  per  acre  foot  of  the 
reservoir  and  a  steam  quality  of  40%  at  the  sandface 
(Bursell  et  al,  1975;  Farouq  Ali,1979  and  Doscher  et 
al,  1979). 

The  first  paper  ever  presented  on  steam  injection  was  done 
by  Ramey,  (1962)  in  which  he  developed  equations  for  the 
estimation  temperature  profile  as  a  function  of  depth  and 
time  for  a  single  phase  flow.  The  modelled  generated  was 
improved  by  (Satter,  1965)by  considering  the  changes  in 
steam  quality,  Overall  heat  transfer  coefficient  and  the  fluid 
properties  which  were  not  accounted  by  Ramey  (1962).  He 
thus  presented  better  equations  that  compared  the  per  cent 
of  heat  loss  for  superheated  steam,  saturated  steam  and 
understated  steam 

(Hoist  &  Rock,  1966)  was  able  to  account  for  the  effect  of 
frictional  loss  and  kinetic  energy  changes  by  dividing  the 
entire  injection  system  into  three-  (a)  flowing  fluid,  (b) 
wellbore,  and  (c)  formation  with  each  part  being  treated 
separately  and  assuming  they  were  interconnected  only  by 
heat  transfer.  This  study  shows  that  steam  quality  can  be 
greatly  affected  by  friction  losses. 

As  an  important  steam  optimization  parameter,  (Willhite, 
1969)  presented  an  iterative  method  for  predicting  the 
overall  heat  transfer  coefficient  by  considering  the  various 
heat  transfer  mechanisms  and  thus  presenting  a  method  for 
the  calculation  of  heat  transfer  coefficient  for  radiation 
through  the  annulus  (hr)  and  the  heat  transfer  coefficient 
for  natural  convection  and  conduction  in  the  annulus  (he). 
Earlougher  (1969)  applied  a  depth-step  technique  similar  to 
Satter’s  for  calculating  heat  losses  and  downhole 
conditions.  He  extended  Satter's  approach  by  including  the 
effects  of  pressure  changes  in  the  injection  tubing  and  the 
effect  of  casing  cement  on  heat  transfer  and  studied  the 
effects  of  various  well  completion  schemes.  He  was  able  to 
demonstrate  the  importance  of  including  the  static  pressure 
term  in  the  pressure  change.  He  also  showed  that  by  using 
insulated  tubing  heat  loss  could  be  reduced  significantly. 
Earlougher  concluded  that  the  bottomhole  properties  of 
steam  are  a  function  of  injection  conditions  and  well 
completion,  and  also  emphasized  that  the  pressure  change 
cannot  be  neglected  in  heat  transfer  calculation  of  steam 
injection.  Pacheco  &Farouq  (1972)  presented  an  analysis  of 
wellbore  heat  losses  and  pressure  drop  for  steam  injection 
assuming  the  steam  to  be  a  perfectly  homogeneous  two- 
phase  flow.  Their  studies  showed  that  an  increase  of 
injection  rate  reduced  heat  loss  and  illustrated  that  frictional 
losses  are  important  in  determining  downhole  steam 
pressure,  quality  and  temperature.  This  study  was  followed 
by  Farouq  Ali  (1981)  developed  a  comprehensive 
mathematical  model  to  simulate  the  vertical  upward  and 
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downward  flow  of  wet  steam  in  a  well.  This  model  is  a 
combination  of  the  Pacecoc£Farouq(1972)  model  and  the 
pressure/flow  regime  correlations  of  Gould  et  al  (1974), 
Chierici  et  al  (1974),  and  Duns  and  Ros  (1961). 

Farouq  Ali  (1986)  used  the  Duns  and  Ros  flow  pattern  map 
to  determine  flow  regime  for  wet  steam  flow.  His  study 
showed  it  is  necessary  to  include  slip  and  flow  regime  in 
calculating  the  pressure  change.  Furthermore,  using  the 
Duns  and  Ros  flow  regime  map,  the  flow  regime  was  found 
to  be  predominantly  in  the  Slug -Froth  flow.  Still  on  flow 
regime  map,  Sylvester  (1984)  has  shown  that  the  Tait elet  al 
(1980)  flow  regime  map  is  superior  to  the  Duns  and  Ros 
map  since  it  predicted  the  annular-mist  flow  at  much  lower 
superficial  gas  velocities  especially  at  higher  pressures. 

As  an  improvement  to  Pacheco  &Farouq  (1972), 
Fontanilla&  Aziz  (1982)  developed  a  mathematical  model 
for  wellbore  heat  loss  that  incorporated  empirical  two- 
phase  flow  correlations  using  Beggs&  Brill  (1973),  Aziz  et 
al  (1972)  and  Yamazaki  &  Yamaguchi  (1979)  correlations. 
Yao  and  Sylvester  (1987)  have  shown  that  the  Beggs  and 
Brill  correlation  is  unsatisfactory  for  vertical  annular-mist 
flow. 

Another  innovative  study  on  steam  injection  was  done  by 
Jianshe^t  al  (2010).  Their  approach  was  able  to  account  for 
the  effect  of  the  reservoir  backpressure  on  the  injection  rate 
and  consequently  the  steam  quality  by  using  a  Nomograph 
developed  for  the  Mukhaizna  Field  as  against  the 
conventional  classical  models  that  neglect  the  inpact  of  the 
reservoir  back  pressure. 

Most  of  the  various  classical  models  above  has  been  used  to 
develop  the  algorithm  used  by  many  steam  injection 
softwares  but  the  most  common  approach  has  been  the 
independent  analysis  of  the  injection  well  and  the  reservoir. 

II.  MODEL  FORMULATIONS 

Considering  a  steam  injection  well  model  as  shown  in 
Figure  (1)  transferring  a  steam-hot  water  mixture  through  a 
control  mass,  AM,  the  general  energy  equation  for  the 
system  at  any  two  unique  conditions  (points)  can  be  written 
as: 


u  _l_  _i_ 

1711  9c'  J 


bni 


=  K 


m2 


+  + 


29cJ  9c  J  29cJ 

In  differential  form,  Equation  (3.1)  becomes 
VmdVm 


—  (2.1) 


dhm  +  — + 

771  9c  J 


29cJ 


M) 


(2.2) 


If  we  include  the  heat  loss  term  and  assuming  no  work  done 


by  or  on  the  steam,  we  have; 

dhm +  ■?-.-+  ^^-dQ=  0  (2.3) 

771  9c  J  29cJ 


Equation  (2.3)  describes  a  general  energy  equation  for  the 
energy  balance  of  the  steam-hot  water  mixture  in  the 
system  For  a  general  concern,  it  is  often  desired  to  express 
the  energy  equation  as  a  gradient  of  the  depth  for  the 
injection  well  optimization.  Hence  we  can  have, 
dhm  _|_  g_  1  _j_  Vm  dVm  _  dQ_  q  (2  4) 

dz  gc  J  gcJ  dz  dz 


z 


Fig.l:  A  Steam  Injection  Well  Model 


Sometimes,  it  is  more  convenient  to  express  the  steam 
mixture  velocity,  Vmin  terms  of  the  superficial  velocity  or 
the  mass  flux  rate  and  specific  volume  of  the  individual 
components  of  the  steam  stream  By  definition, 

Kn  ~  Vsl  +  Vsg  ~  Glvl  +  Ggvg  (2-5) 

dVm  =  Gi  dvt  +  Gg  dvg  (2.6) 

Where; 


Vsl  =  liquid  (hot  water)  superficial  velocity,  ft/hr 
Vsl  =  gas  (vapour)  superficial  velocity,  ft/hr 
Gt  =  liquid  mass  flux  rate,  lb/hr-ft2 
Gg  =  gas  mass  flux  rate,  lb/hr-ft2 
vz  ^liquid  specific  volume,  ft3/lb 
vg  =  gas  specific  volume,  ft3/lb 

Putting  (2.4)  and  (2.5)  into  (2.6),  we  have  the  following 
result; 


Hhn+  JLI+  —  [v;G;2  —  4-  vlGlGgd-^+  vgGgGl^  + 

dz  dr  I  drjy-  L  L  dz  L  L  9  dz  9  9  L  dz 

(2.7) 


9c 

Gq  ““4  —  — 

9  dz  J  dz 


Some  of  the  essential  properties  of  the  steam  that  is  of 
primary  interest  to  this  study  are  the  mixture  enthalpy,  gas 
specific  volume  and  liquid  specific  volume  defined  as; 

h™  =  fiX,  P),  hence,  ^^-can  be  evaluated  as  follows; 

"  dz 

dhm  dhm  dX  ^  dhm  dP  ^ 

dz  dX  dz  dP  dz 
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More  so, 


dva 


dvg  dP 


(2.8b) 


(2.8c) 


dz  dP  dz 
dvi  _  dvi  dP 
dz  dP  dz 

Substituting  Equations  (3.8a-c),  and  solving  for  steam 
quality  gradient,  we  can  have  the  following  result; 

dX  _  1  r  dhm  dP  g  1  1  dP  f  2  dvt 

dz  L  dp  dz  9c  J  9cJ  dz\  1  1  dP 


dX 
dv . 


v,G,G  —  +  vgGgGl—+  vGG2g— )-— 1 

1  1  9  dp  9  9  I  QP  G  9  gp  J  dz  j 


(2.9) 


avg 

dP  a  -  y  *  dP  u  9  dp  . 

The  Equation  (2.9)  above  is  a  first  order  differential 
equation  that  can  be  used  to  estimate  the  steam  quality 
gradient  analytically. 

dP 

To  estimate  the  pressure  drop  term, — ,  we  establish  the 

dz 

momentum  balance  equation  in  terms  of  the  mechanical 
energy  balance  of  the  system  This  is  also  sometimes 
conventionally  express ed  as  the  pressure  drop  equation  for 
fluid  flow  through  a  pipe  section.  Hence  we  can  write  that; 

(2.10a) 


dP 


- —  dz  +  VjD^nL  +  dWs  +  dWf  =  0 

Pm  9c  9c 


Where: 

dP=  total  Pressure  differential,  lb/ft2  (Psf) 
dWs  =  Work  done  by  or  on  the  fluid,  lb -ft/lb 
dWf  =  Frictional  work,  lb  -ft/lb 
For  dP'm  Psi,  we  can  rewrite  (3.10a)  as; 

144dP  _  3_dz  +  YmMm  +  ^  +  dW '  =  0  (2.10b) 

Pm  9c  9c 

Since  there  is  no  Work  done  on  or  by  the  steam  stream 

dWs 

(expansion  or  compression  of  the  fluid  (steam)),  — »  0  such 
that  Equation  (3.10b)  becomes, 


144  dP 


Pm9  ,  Vm  dVm  ,  _ 

- 1 - b  Pr 

9c  9cdz 


dWf 


=0 


(2.11) 

dz  gc  gcdz  dz 

As  usual  by  solving  for  the  total  pressure  gradient,  we  can 
establish  this  equation  below; 


d_P 

dz 


144+ 


Pm\ 

9c 


h  c? 


Pm9__fdP\ 

9  c  ''dz'  frjctionlosses _ 

dv,  dv  q  dv,  ndvq 1 

Jf+  nGlGgjf+  ''gGgG,—+  vcGljfl 


(2.12) 


In  this  study,  the  steam  properties  and  the  injection  well 
conditions  will  be  generated  using  PROSPER  and  a 
sensitivity  test  will  be  run  using  critical  parameters  as 
presented  in  the  next  chapter. 

1.  The  Reservoir  Back -Pres sure  Effect 
It  is  a  common  experience  that  the  reservoir  pressure  causes 
a  significant  constraint  during  steam  injection.  Therefore,  a 
total  neglect  of  this  phenomenon  will  limit  the  accuracy  of 
the  predicted  steam  properties.  The  back-pressure 
phenomenon  can  be  modeled  using  the  figure  below 
(Figure  2). 

By  capillary  effect,  the  in-situ  fluid  tends  to  rise  through  the 
vertical  column  of  the  injector.  For  this  set  up,  if  we  assume 
that  the  steam  generator  discharge  pressure  remains 


unchanged  at  the  wellhead,  the  total  pressure  of  the 
injection  well  system  at  a  constant  injection  rate  can  be 
established  thus; 

Pinj  +  ( Pres  ~  0.052 Podh)  =  P  (2.13) 

Where; 

P=  Actual  Downhole  Steam  Pressure  (which  is  equivalent 

to  the  total  pressure  of  the  system  under  injection 

conditions),  Psi 

p0  =  Oil  density,  ppg 

Pinj  =  Injection  Pressure  at  depth  dz,  psi 


Fig.2 :  A  well  Model  for  Reservoir  Back  Pressure  Effect 


It  is  worth  noting  that  this  pressure  ( Pinj- )  at  the  sandface 
(dz=Z)  can  be  greater  or  less  than  the  injection  constraint  at 
the  wellhead  (dz^O)  depending  on  the  dominating  factor 
during  the  steam  transfer  to  the  downhole.  If  gravity 
dominates  flow,  Piri >  Piri but  when  frictional 

drag  dominates  the  steam  flow,  \Pirii  ]  <  \Pirii  ]  .  It 

only  becomes  unchanged  when  both  the  gravity  losses  and 
the  friction  losses  have  approximate  equal  inpacts  on  the 
injection  system 

As  ‘dh  ’  approaches  zero  TVD,  the  influence  of  the 
reservoir  pressure  due  to  the  capillary  column  of  the  rising 
becomes  less  significant. 

In  practice,  the  injection  well  is  totally  filled  with  the  steam 
column  such  that  ^approaches  Z.  This  is  certainly  the  case 
for  heavy  oil  wells  which  do  not  readily  flow  by  natural 
effect  and  also  as  a  result  of  external  constraint  of  the 
injection  pressure.  Therefore,  we  can  rewrite  Equation 
(2.13)  as, 
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p.  .  +  p  =  p 

1  in  j  1  1res  1 


(2.14) 


In  other  words; 

Pinj  (l  +  ^s)  =  Pin]  (1  +  iPP)  =  P  (2.15) 

Where  ipPis  defined  as  the  pressure  ratio  of  the  reservoir  to 
the  Steam  injection  pres  sure. The  term  (l  +  i/jp )  is  defined 
as  the  reservoir  back  pressure  (RBP)  correction  factor 
denoted  in  this  study  as  ‘%P\  Therefore,  the  corrected 
pressure  of  the  steam  system  can  be  express  ed  as, 

P  =  PlnJ-$p  (2.16) 


Fig. 3:  Algorithm  for  Reservoir  Back  Pressure  Analysis 


The  implication  of  the  above  Equation  (2.16)  is  that  a  less 
accurate  steam  quality  (overestimation)  will  be  made  by 
ignoring  the  reservoir  back  pressure  especially  for 
relatively  high  pressure  reservoirs  as  this  tends  to  constrain 
the  injection  pressure  and  consequently  the  injection  rate  of 
the  steam  For  convenience  sake,  this  equation  will  be  used 
to  establish  a  table  of  values  of  P  for  different  pressure 
ranges  as  presented  in  the  appendix  section.  The  effect  of 
reservoir  back  pressure  on  the  steam  properties  and  the 
injection  well  system  was  analyzed  numerically  using 
PROSPER  simulator  and  the  result  presented  in  the  next 
chapter.  To  achieve  this,  the  reservoir  pressure  was  set  at 
the  corrected  pressure  based  on  Equation  (3.16)  using  the 
following  algorithm 

2.  The  Wellbore  Heat  Loss  Calculations 

In  the  literature  review  of  this  study  as  presented  in  the 
preceding  chapter,  it  is  clearly  identified  that  the  most 
influencing  factor  for  the  optimization  of  steam  injection  is 
the  choice  of  completion.  This  is  because  the  completion 
design  directly  affects  the  heat  losses  that  occur  in  the 
injection  well.  A  typical  model  for  this  is  given  in  the 
figure  below  (Figure  4) 

I 


Fig.4:  A  Steam  Injector  Heat  Loss  Schematic  Based  on 
Tubing -Inserted-in-Tubing  Model 


Using  the  figure  above,  the  heat  transfer  mechanism  by 
conduction,  convention  and  radiation  can  be  modeled  based 
on  the  following  considerations/as  sumptions. 

•  Heat  transfer  in  the  injection  well  system  is  at 
pseudo-steady  state 

•  The  tubing  and  the  casing  are  symmetrically 
vertically  placed 

•  There  is  no  annular  refluxing 

•  The  overall  heat  transfer  coefficient  is  independent 
on  depth 
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If  we  proceed  with  these  assumptions,  the  following  heat 
transfer  summary  can  be  written  for  the  Model  above  in 
Figure  (4). 

•  The  heat  transfer  to  the  inner  tubing  wall  due  to 
steam  motion  is  by  convention 

•  The  heat  transfer  between  the  inner  tubing  wall 
and  the  outer  tubing  wall  is  by  conduction 

•  The  heat  transfer  between  the  outer  tubing  wall 
and  the  solid  insulation  (for  insulated  tubing)  is  by 
conduction 

•  The  heat  transfer  between  the  insulator  and  the 
annular  space  is  by  the  annular  fluid  convention  and 
conduction  (if  there  is  annular  fluid)  and  the  insulator 
radiation 


Overalltemperature  difference  =  —  Te  (3.21) 
Hence,  the  final  equation  for  the  overall  heat  transfer  of  the 
steam  injection  system  becomes; 

J  =  2nrt0Ut0(Tf-Te )  (3.22) 

This  Equation  can  be  equated  to  the  Ramey’s  premier  heat 
loss  per  unit  length  of  injection  path  given  as; 

2nr‘-iST’ -T*)  ■  R23> 

Where; 

Rt0  =  overall  thermal  resistance 

During  steam  injection,  the  compensation  for  heat  losses  is 
basically  by  reducing  Ut0  as  low  as  possible  by  both  tubing 
and  annular  insulation.  If  we  therefore  solve  for  Ut0  using 
Equations  (3.20a-f)  through  (3.23)  it  will  give; 


•  The  heat  transfer  between  the  annulus  and  the 

inner  casing  wall  is  by  conduction 


•  The  heat  transfer  between  the  inner  casing  wall 
and  the  outer  casing  wall  is  by  conduction 

•  The  heat  transfer  between  the  outer  casing  wall 
and  the  cement  bond  is  by  conduction 

•  The  heat  transfer  between  the  cement  bond  and  the 
adjacent  formation  is  by  conduction 

From  Equations  (2.4),  (2.5)  and  (2.6),  the  1-D  heat  transfer 
Model  through  a  hollow  cylinder  can  be  written  thus; 


dq 

dz 


2nk 


(7)  -  T0){f  orheatconduction}  (2.17) 


=  ZnrihtTi  —T0){forh  eat  convent  ion]  (2.18) 
—  =  Znr^iJi  -  T0){forheatradiation }  (2.19) 


From  the  heat  transfer  summary,  we  can  solve  for 
temperature  differences  as  follows  so  as  to  define  the 
overall  heat  transfer  coefficient 


Uto  = 


kt  kins 


rto^  rtolnl«. 


(3.24) 


( hc+hr ) 


For  a  steel  casing  and  tubing  (or  any  other  high  conductive 
metals), 

kt  =  kc  »  kins  ,  kcem,  hcandhr  .  Therefore,  Equation 
(3.24)  can  be  simplified  to  the  following; 

Uto  = 


rto  In-J 


£fi__ 


(hc+hr) 


.2^ 


(3. 24.  a) 

If  no  insulation  of  the  tubing 

(3.24.b) 


Ut0  = 


rtolnFA 


tinf 


(hc+hr) 


For  no  annular  insulation  and  tubing  insulation,  Ut0  = 

7^?  (3‘24c) 


The  terms,  he  and  hr  can  be  evaluated  independently  based 
on  empirical  correlations/  charts  for  the  different  choice  of 
insulation  material  as  presented  in  Fidan  (2011)  studies. 
Also  using  PROSPER  wellbore  simulator  by  selecting 
Enthalpy  Balancing  Model,  the  overall  heat  transfer  of  the 
injection  well  can  be  calculated  if  there  is  adequate  data  for 
formation  lithology  and  the  completion  status . 

For  the  purpose  of  this  study,  in  order  to  numerically 
analyze  the  inpact  of  injection  well  completion,  a 
sensitivity  test  was  run  for  known  values  of  overall  heat 
transfer  coefficient. 


Using  the  Equations  (3.26a-f),  the  overall  temperature 
difference  becomes 
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3.  Study  Simulation  Methodology 


Fig.  5:  Simulation  Flow  Chart 


III.  MODEL  RESULTS  AND  DISCUSSIONS 

The  analysis  of  steam  injection  starts  with  the  specification 
of  the  input  steam  data  and  then  the  injection  well 
configuration.  A  special  consideration  was  made  for 
injecting  at  the  same  tempeerature  and  different  overall  heat 
transfer  coefficient  and  injecting  at  the  same  overall  heat 
transfer  coefficient  and  different  temperatures . 

Based  on  the  specified  data  of  the  input  steam  and  the 
injection  well  configuration  used  in  this  study,  the  steam 
data  generated  at  a  pressure  of  llOOpsiwas  presented  in  Fig 
6  and  Fig  7. 


(a)  Injection  Temp  =  700°F  (b)  Injection  Temp =600° F 

Figure  6:  Steam  Quality  Gradient  at  Different  Injection 
Temperature  and  Uto 


(a)  Injection  Temp=700°F  (b)  Injection  Temp  =  700°F 
Fig. 7:  Steam  Temperature  Gradient  at  Different  Injection 
Temperature  and  Uto 


From  the  figuers  above,  it  can  be  noted  that  though 
temperature  of  700°F  gave  a  good  SQ  but  in  as  much  as  a 
sufficient  steam  quality  can  be  generated,  it  will  be  more 


www.iiaers.com 


Page  |  116 


International  Journal  of  Advanced  Engineering  Research  and  Science  (IJAERS)  [Vol-5,  Issue-11,  Nov- 2018] 

httos:/ Zdx.doLora/1 0.221 61/iiaers.5.1 1.16  ISSN:  2349-6495(P)  /  2456-1908(0) 


economical  to  consider  injecting  at  600°F.  The  differences 
in  the  Overall  heat  transfer  coefficeint  was  used  to  analyse 
the  possible  necessity  of  insulating  the  wellbore  or 
inproving  the  insulation  efficiency.  Using  this  as  a  design 
guide  shows  that  an  overall  heat  transfer  coefficient  of 
8Btu/hr/Ft2/°F  can  be  adequate  for  this  operation  since  it 
gave  a  sufficient  SQ  used  to  specify  the  Injection  well 
control  of  the  ECLIPSE  programme. 


(a)  Injection  Temp  =  700°F  (b)  Injection 

Temp=600°F 

Fig. 8:  Steam  Pressure  and  Enthalpy  Gradients 

1.  Reservoir  Back  Pressure  (RBP)  Sensitivity 
Analysis 

A  major  limitation  to  classical  wellbore  SQ  analysis  is  total 
neglect  of  the  effect  of  reservoir  back  pressure  on  the  steam 
injection  rate  and  consequently  the  SQ.  The  result  in  Table 
(1)  is  a  sensitivity  study  performed  with  a  5000ft  steam 
injector  where  the  impact  of  reservoir  backpressure  is  more 
pronounced.  Hence,  it  was  observed  that  the  reservoir  back 
pressure  causes  a  significant  constraint  to  the  injection  rate 
with  a  resultant  significant  SQ  drop.  Therefore,  an  optimal 
performance  of  steam  injection  design  will  require  a 
negligible  reservoir  pressure.  The  phenomenon  of  RBP  is 
also  graphically  illustrated  in  Figure  (9)  using  a  simple 
MATLAB  multiple  plot  tool. 


Table. 1:  SQ  Vs  Depth  at  Different  Reservoir  Pressure  (RBP 
Sensitivity  Analysis) 


DEPTH 

(ft) 

SQ  @ 
fp=  1 

SQ  @ 

=  1. 1429 

SQ  @ 

=  1.2857 

SQ  @ 

%p 

=  1.4286 

SQ 

@ 

%p 

=  1.  57 

0 

100 

100 

100 

100 

100 

250 

100 

100 

100 

100 

100 

500 

100 

100 

100 

100 

100 

750 

100 

100 

100 

100 

100 

1000 

100 

100 

100 

100 

100 

1250 

100 

100 

100 

98.22 

94.96 

1500 

100 

99.3 

97.07 

94.14 

90.22 

1750 

98.22 

96.14 

93.51 

90.07 

85.51 

2000 

95.37 

92.98 

89.96 

86.02 

80.8 

2233.3 

92.75 

90.08 

86.69 

82.3 

76.42 

2466.7 

90.22 

87.25 

83.5 

78.65 

12.1 

2700 

87.75 

84.5 

80.4 

75.05 

67.82 

2933.3 

85.35 

81.81 

77.36 

71.5 

63.6 

3166.6 

83.02 

79.22 

74.37 

67.99 

59.41 

3400 

80.75 

76.68 

71.41 

64.51 

55.26 

3633.3 

78.57 

74.18 

68.5 

61.09 

51.14 

3866.6 

76.43 

71.72 

65.62 

57.7 

47.05 

4100 

74.33 

69.3 

62.78 

54.34 

42.99 

4333.3 

72.28 

66.91 

59.98 

51.01 

38.96 

4566.6 

70.26 

64.56 

57.21 

47.71 

34.94 

4800 

68.27 

62.24 

54.48 

44.4 

30.94 

4900 

65.91 

59.5 

51.28 

40.66 

26.38 

5000 

NOTE 

63.57 

56.79 

48.11 

36.9 

21.8 

Liqui 

d 

Well 

head 

Liquid 

Rate 

Liquid 

Rate 

Liquid 

Rate 

Liquid 

Rate 
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@RBP 
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ts 

=1400Psi 
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=0Psig 
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400Psig 
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P= 

8  OOP 
sig 

1217. 

Uto=5Bt 
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1412.4 

8 

u/hr/ft2/F 
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Fig.9:  Reservoir  Back  Pressure  Sensitivity 


2.  Reservoir  Grid  Definition 

Reservoir  X  is  a  9x5x4  reservoir  with  active  cells  specified 
in  X-Y  plane.  There  are  Three  wells-  an  injector  at  (1,1) 
and  two  producers  at  (5,5)  and  (9,1)  respectively  with 
different  transmissibities  at  each  Z-layer.  An  injection 
design  was  specified  for  10  years  injection  and  a  timestep 
of  365  days  for  selected  for  analysis. 


Fig.10:  Reservoir  Grid  Layerout  in  3-D 


- 6.T3E  (BMLEnaPCflr.BB) 


Fig.  11 :  Field  Presure  History  for  the  lOyrs  injection 
Period 


3.  Field  Pressure  Response 

The  Figure  11  above  displays  the  pressure  response  of  the 
steam  injection  from  an  initial  datum  depth  pressure  of 
75psi  to  a  stablized  pressure  of  about  1090psi  which  shows 
that  steam  does  not  just  add  thermal  energy  to  the  reservoir, 
hence,  it  also  provides  a  water  drive  to  the  reservoir  thereby 
aiding  the  area  sweep  efficiency. 

4.  Field  Heat  Loss  Rate  (FHLR)  and  Heat  Loss 
Total  (FHLT) 

When  steam  is  finally  introduced  into  the  reservoir  at  a 
particular  SQ,  the  efficiency  of  the  heavy  oil  displacement 
will  be  characterized  by  the  rate  of  heat  loss  throughout  the 
injection -period. 

Figure  12 -(a):  Field  Heat  Loss  Rate 
Figure  12-(b):  Field  Heat  Loss  Total. 

Figure  1(a)  shows  the  simulated  field  rate  of  heat  loss 
behavoiur.  It  can  be  clearly  seen  that  the  FHLR  increses  at 
the  initial  period  of  injection  of  about  2years  (700days)  of 
injection  before  a  decline  in  the  FHLR.  This  is  because  at 
earlier  stage  of  injection,  the  rate  of  heat  conduction  to  the 
overlaying  formation  is  high  and  it  is  propagated  similarity 
as  a  pressure  transient.  Hence,  more  heat  is  needed  to  warm 
up  the  reservoir  than  to  keep  it  at  a  stablized  temperature. 
This  is  more  clearly  illustrated  in  Figure  12(b)  which  shows 
that  at  a  later  stage  of  the  injection,  FHLT  becomes  near 
constant.  The  resulting  field  temperature  profile  during  the 
periods  of  injection  is  shown  in  Figures  13(a)  to  13(c). 
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Fig, 13(a):  Field  Temperature  Profile  Before  Steam 
Injection  ( 125°F) 

Fig  13(b)  shows  that  after  lyr  of  injection  both  PROD1  and 
PROD2  has  not  been  adequately  heated  up  while  after  2yrs, 
PROD1  has  been  well  affected  by  the  steam  injection  while 
PROD2  remains  slighly  affected  as  shown  in  Fig  13(c). 
This  is  because  PROD1  is  closer  to  the  INJ1  since  heat 
distribution  is  both  time  and  space  dependent. 


k = ■  a  jr  m 


frWilM  h 


Fig.13:  (b)  Field  Temp.  Profile  after  1  yr  (c)  Field  Temp. 
Profile  after  2  yrs 

5.  Field  WaterProduction  and  Field  Oil  Production 
Analysis 

A  comon  experience  with  steam  injection  project  is  the 
increased  water  cut  as  the  water-oil  front  changes  with  time 
and  space  and  therefore  tends  to  the  producers.  As  a  result 
of  excessive  turbulence  in  the  reservoir  caused  by  the 
injection,  the  condensed  steam  in  the  hot  water  zone  are 
produced  along  with  the  heated  oil  bank.  Henec,  water 
production/watercut  increases  with  time.  A  good  advantage 
of  this  process  is  that  it  results  to  a  secondary  water  drive 
mechanism  that  effects  the  sweeping  of  the  heated  oil  bank. 
Figures  14(a)  to  14(c)  shows  the  field  water  production 
profile  and  the  individual  contibutions  of  the  two 
producers -PROD1  and  PROD2. 
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15(a)  to  15(c).  It  is  observed  that  PROD1  responds  faster  to 
water  cut  than  PROD2  which  responds  slower.  This  can  be 
attribited  to  the  differences  in  the  location  of  the  Producers 
with  respect  to  the  Injectors.  The  steam  front  during 
injection  reaches  PROD1  earleir  than  PROD2.  The 
implication  of  this  as  illustrated  in  the  simulation  results  is 
that  subsequently  at  a  much  later  period,  the  PRODlwillbe 
dominated  by  hot  water  zone. 


A  comparative  Analysis  of  Field  Oil  Production  Rate -Field 
Water  Production  Rate  and  Field  Oil  Production  Total- 
Field  Water  Production  Total  are  given  in  Figures  16(a)  and 
16(b).  The  plots  shows  that  as  water  production  rate 
increases,  oil  production  Rate  decreases.  The  Figure  16(b) 
precisely  displays  a  saturation  rate  growth  model  and  a  near 
linear  relationship  for  the  oil  production  and  water 
productions  Total(s)  of  the  injection  period  respectively. 


( ii)  Field  Water  Production  Rate 
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producers.  The  counter  result  of  this  process  is  the  increase 
in  water  saturation.  Figure  17(a)  shows  the  oil  saturation 
history  of  the  reservoir  for  lyear  (365  days),  5years  (1825 
days)  and  3years  (3650days)  of  injection  respectively. 

In  the  literature  review,  it  was  clearly  stated  that  steam 
injection  activates  solution  gas  drive  in  the  reservoir  system 
and  also  improves  the  quality  of  the  produced  oil  by 
thermal  distillation.  Hence,  steam  injection  increases  the 
gas  saturation  as  the  input  thermal  energy  of  the  steam 
helps  to  liberate  the  light  gas  molecules  in  solution.  The 
economic  benefit  of  this  is  that  recovery  by  thermal 
methods  upgrades  the  API  gravity  of  the  in -situ  oil.  Figure 
(17)  shows  the  oil  saturation  history  while  the 
corresponding  gas  saturation  history  of  the  reservoir  is 
shown  in  Figure  (18). 


6.  Effect  of  Steam  Injection  on  Oil  and  Gas 
Saturations 

Steam  injection  like  water  injection  causes  the  in -situ  oil 
saturation  to  vary  with  both  time  and  space  as  the  steam 
front  propagates  along  the  reservoir  grid  cells  towards  the 


365  Days  1825  Days  3650  Days 


(a)  (b)  (c) 

Fig. 17:  Field  Oil  Saturation  History 
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FigJ8:  Field  Gas  Saturation  History 
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IV.  CONCLUSION 

This  study  presents  an  integrated  steam  injection  Model 
needed  for  the  complete  optimization  of  the  steam  injection 
project.  The  study  frame  work  incorporates  a  well  computer 
simulator  known  as  PROSPER  (An  IPM  Suit  Product)  and 
a  reservoir  simulator  known  as  ECLIPSE  (Eclispse300)  for 
the  analysis.  Thus,  the  optimal  use  of  the  injected  fluid  was 
archived  since  a  semi-iterative  procedure  can  be  used  for 
data  convergence  as  shown  in  simulation  flow  chart. 

The  steam  properties  predicted  with  this  study  shows  that 
higher  SQ  will  be  generated  at  higher  injection  Pressure 
and  Temperature  and  at  lower  overall  heat  transfer 
coefficient. 

By  adequately  reviewing,  formulating  and  analysing  the 
basic  analytics  behind  the  steam  injection  mechanism,  the 
impact  of  reservoir  back  pressure  on  the  SQ  and 
consequently  the  injection  rate  was  numerically  quantified 
using  a  5000ft  case  study  well  at  a  constant  injection  well 
head  constraint.  The  results  showed  that  a  most  efficient 
steam  injection  will  be  achieved  at  the  lowest  reservoir 
initial  pressure. 

The  various  sensitivity  results  generated  with  ECLIPSE 
showed  that  production  well  spacing  relative  to  the  injector 
have  a  critical  effect  on  the  area  sweep  efficiency  of  the 
injected  steam  especially  in  the  early  stage  of  injection  and 
hence  optimal  well  spacing  should  also  form  the  basic 
development  criteria  for  a  given  field. 


1 .  Recommendati  ons 

The  efficiency  of  steam  injection  is  basically  dependent  on 
the  effective  SQ  generated.  To  generate  a  high  SQ 
especially  in  relatively  deep  wells,  the  only  option  has  been 
to  effectively  thermally  insulate  the  tubing  and  the  annulus 
which  will  make  the  steam  injection  project  a  non -attractive 
venture  because  of  the  increased  cost  of  insulation.  To 
compensate  for  this,  a  consideration  can  be  made  for 
sourcing  the  insulation  materials  locally  using  the  materials 
in  the  table  below: 

Table.2:  The  Thermal  Conductivities  of  Some  Locally 


Sourced  Insulation  Materials 


Material 

Thermal 
Conductivity 
(W/m  k) 

Thermal 
Conductivity 
(Btu/Ft  °F) 

Shredded 

0.17 

0.1115 

Asbestos 

Sheets 

Dry  Ash 

0.12 

0.078744 

Cork,  Felt 

0.05-0.10 

0.0328  - 

0.0.06562 

Freon 

0.0083 

0.0054465 

Moreover,  the  difficulty  involved  in  integrating  two  unique 
simulators  made  the  study  more  tedious.  Based  on  this,  the 
study  only  forms  a  framework  of  a  future  integrated  steam 
injection  simulator  which  can  be  more  efficient.  Hence,  this 
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research  work  welcomes  any  idea  on  the  development  of  a 
fully  integrated  steam  injection  model/simulator. 

As  part  of  the  above  recommendation,  such  proposed 
simulator  should  be  able  to  auto -calculate  the  overall  heat 
transfer  coefficient  based  on  the  defined  completion  status. 
This  is  necessary  to  avoid  unnecessary  switching  between 
the  rigorous  analytics  involved  in  the  estimation  of  the 
overall  heat  transfer  coefficient  that  often  requires  a  third 
party  simulator. 

2.  Contributions  to  Knowledge 

This  study  provides  a  frame  work  for  the  complete 
optimization  of  steam  injection  design  and  as  such 
constitutes  a  novel  approach  to  steam  injection  modelling. 
Also  very  important  to  mention  is  the  novel  approach  to  the 
analysis  of  Reservoir  Back  Pressure  impact  on  Steam 
Quality  predicted  during  the  steam  injection  processes. 
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